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TABLE  1.  Statistics  and  their  control  limits  in  modified  X  and  T 
for  monitoring  the  means  of  a  multidimensional  production  process. 


FIGURE  1. 


Geometric  distributions  of  run  lengths  for  the  X  and  T' 
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monitoring  at  level  a  =  0.05  for  the  case  p*2,  n*6  ,  with  6=3.0  and  6=6.75. 
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0.  Abstract.  Hotelling's  T  and  related  charts  are  studied  in  moni- 

2 

toring  the  means  of  a  multidimensional  production  process.  The  T  chart 
is  shown  to  be  optimal  in  a  class  of  procedures  in  that  it  signals  more 
quickly  than  other  procedures  in  the  class  when  the  process  is  not  in  con¬ 
trol.  Nonstandard  properties  of  the  distributions  of  run  lengths  of  these 
charts  are  studied  when  (i)  certain  parameters  are  estimated  in  a  base 
period  and  modified  procedures  are  followed  using  these  estimates,  (ii) 
the  process  is  a  drifting  process,  and  (iii)  the  assumption  of  independent 
Gaussian  vector  observations  is  replaced  by  the  assumption  that  observa¬ 
tions  are  generated  from  a  spherical  process.  For  these  cases  stochastic 
bounds  on  the  actual  run-length  distributions  are  given  in  terms  of  geo¬ 
metric  distributions,  and  certain  monotone  properties  of  run  lengths  are 
established  under  drifting. 
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1.  Introduction.  Various  control  charts  are  available  for  monitoring 

the  output  of  a  production  process.  Shewhart's  (1926)  X  chart  is  standard 

for  monitoring  the  mean  level  of  a  single  quality  characteristic,  and 

2  2 

Hotelling's  (1947)  X  and  T  charts  are  used  for  monitoring  the  means  of 

several  characteristics.  In  these  charts  successive  values  of  a  statistic 

are  plotted  against  time,  and  a  chart  signals  at  level  a  that  the  process 

is  not  in  control  when  the  statistic  exceeds  a  control  limit  c  .  Standard 

a 

assumptions  are  that  independent  random  samples  of  n  observations  each  are 

taken  in  succession  from  a  p-dimencional  Gaussian  process  having  the  mean 

vector  u  and  the  dispersion  matrix  E,  and  is  the  value  of  u  when  the 

2 

process  is  in  control.  For  the  X  chart  the  statistics  are 

xi a  n(Ti-y0>^"1(fi-y0)’  1  =  (l.D 

where  is  the  vector  of  means  from  the  ith  sample,  and  the  control  limit 
2 

is  c^  =  xa(p)»  the  100O-~a)  percentile  of  the  central  chi-squared  distri- 

2 

bution  having  p  degrees  of  freedom.  For  the  T  chart  the  statistics  are 

Ti  "  1  =  d. 2) 

where  is  the  sample  dispersion  matrix  for  use  when  Z  is  unknown  and 

2 

n  >  p,  and  the  control  limit  is  c  *  Ta(p,n-1),  the  lOO(l-a)  percentile  of 
Hotelling's  (1931)  distribution  having  the  parameters  p  and  n-1.  These 
procedures  were  developed  further  and  applied  by  Jackson  (1956,1959), 
Jackson  and  Morris  (1957),  Jackson  and  Bradley  (1961a, b),  Ghare  and  Torg- 
ersen  (1968),  and  Jackson  and  Mudholkar  (1979),  including  diagnostic  pro¬ 
cedures  for  assigning  probable  causes  when  a  chart  signals  that  the  moni¬ 
tored  process  is  not  in  control. 
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Our  purpose  here  is  to  study  the  X  and  T  charts  and  some  modifica¬ 
tions  of  these  under  standard  and  nonstandard  assumptions.  To  be  effect¬ 
ive  a  chart  should  signal  infrequently  when  the  process  is  in  control  and 
more  frequently  as  the  process  shifts  from  control.  These  properties  are 
embodied  in  the  distribution  of  the  run  length  of  a  chart,  i.e.  the  number 
of  successive  samples  taken  before  the  chart  signals  that  the  process  is 
not  in  control.  Run-length  distributions  are  useful  in  comparing  monitor¬ 
ing  procedures,  and  their  moments  enter  naturally  into  cost  models  for 
designing  economic  inspection  policies.  In  Shewhart  charts  successive 
statistics  typically  are  independent;  if  in  addition  the  process  is  sta¬ 
tionary,  then  the  distribution  of  run  lengths  is  geometric  with  parameter 

equal  to  the  probability  of  exceeding  the  control  limit  on  any  sampling 

2  2 

occasion.  Under  standard  assumptions  the  run  lengths  of  the  X  and  T 

charts  are  geometric,  assuring  that  these  charts  eventually  will  signal 

with  unit  probability.  Following  the  supporting  materials  of  Section  2, 

2 

we  develop  further  properties  of  the  standard  T  charts  in  Section  3  in¬ 
cluding  optimal  properties. 

2  2 

Limitations  are  inherent  in  the  standard  X  and  T  charts,  and  these 

inhibit  use  of  the  charts.  Process  dispersion  parameters  frequently  are 

2 

unknown  for  use  in  the  X  charts,  and  the  control  value  uQ  may  be  unknown 

2 

for  use  in  either  chart.  Moreover,  the  T  chart  requires  that  n  >  p  in 
each  sample  to  ensure  the  invertibility  of  S^,  whereas  in  practice  small 
but  frequent  samples  routinely  are  taken.  In  the  case  of  a  single  char¬ 
acteristic  a  common  practice  is  to  estimate  such  parameters  in  a  base 
period  of  sufficient  length  to  give  reliable  estimates,  and  to  use  these 
estimates  in  lieu  of  the  unknown  parameters.  In  the  case  of  several 
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2  2 

characteristics  this  leads  to  modified  X  and  T  charts  as  studied  in  Sec¬ 
tion  4.  One  consequence  is  that  run  lengths  of  these  charts  are  not  geo¬ 
metric  owing  to  dependencies  among  the  successive  statistics. 

2  2 

Nonstandard  properties  of  the  X  and  T  charts  are  examined  further  in 
Sections  5  and  6.  Section  5  deals  with  drifting  processes;  their  run  lengths 
are  not  geometric  because  successive  statistics  are  not  distributed  identi¬ 
cally.  Nonetheless,  a  basic  monotonicity  property  of  the  run  lengths  is 
established  and  used  to  bound  the  actual  run-length  distributions  in  terms 
of  geometric  distributions.  In  Section  6  the  Gaussian  assumption  is  weak¬ 
ened  to  the  case  that  observations  are  generated  from  a  spherical  process, 
and  these  charts  are  studied  under  the  weakened  assumptions. 

2.  Preliminaries.  Let  1/  be  a  finite-dimensional  linear  space;  ex¬ 
amples  are  the  n-dimensional  Euclidean  space  Rn  and  the  space  of  real 

(nxm)  matrices.  Of  two  random  variables  U  and  V,  U  is  said  to  be  stochas¬ 
tically  larger  than  V  if,  for  every  t,  P(U  >  t)  >_P(V>t).  This  in  turn 
implies  that  E(U)  >_  E(V)  whenever  the  expected  values  are  defined.  Let 
{G(t;a);  a  e  (0,1)}  be  cumulative  distribution  functions  (cdf's)  of  geo¬ 
metric  distributions  having  argument  t  and  parameter  a;  then  the  family 

(G(t;a);  a  e  (0,1)}  is  stochastically  decreasing  in  a.  If  u(*)  and  v(*) 

N 

are  two  probability  measures  on  R  ,  u(*)  is  said  to  be  more  peaked  about 
N 

9  e  R  than  v(«)  in  the  sense  of  Sherman  (1955)  if,  for  every  compact  con- 
N 

vex  set  E  c  R  symmetric  about  9  under  reflection,  u(E)  >_v(E). 

A  basic  inequality  is  the  following;  a  standard  proof  uses  convexity 
and  Jensen’s  inequality. 

LEMMA  1.  Let  (X^.X^ , . . . ,Xt}  be  mutually  independent  random  elements 
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with  values  in  l/g  *  l/C  such  that  {X^  ,X£ , .  •  •  ,Xt }  are  distributed  indentically. 
For  any  function  4>:  ^  +  ^  and  any  constant  c, 

t 

P[<KX1,X0)<c,...,<KXc,X0)<c]  1  JjP  [♦(Xi>X(J)<c]. 

The  available  quality  control  procedures  for  variables  are  based  al¬ 
most  exclusively  on  Gaussian  assumptions.  To  study  the  effects  of  depart¬ 
ures  from  normality  we  consider  distributions  having  the  weaker  property 
of  spherical  symmetry,  including  such  heavy-tailed  distributions  as  the 

N 

spherical  Cauchy  law  and  other  spherical  stable  distributions.  If  Y  e  R1, 

is  a  random  element  having  the  distribution  L( Y) ,  then  L( Y)  is  said  to  be 

spherically  symmetric  if,  for  any  (NxN)  orthogonal  matrix  Q,  L(Q Y)  =  L( Y) . 

Properties  of  these  and  related  distributions  were  considered  by  Kelker 

(1970).  Let  be  the  class  of  X-dimensional  spherical  distributions 

symmetric  about  9  e  RN,  i.e.  if  L(Y)  s  SN(9,IN>,  then  L( Y-6)  is  spherical. 

If  M  is  a  linear  subspace  of  R  ,  a  function  d> :  R  -*■  R  is  said  to  be  loca- 

N 

tion- invariant  with  respect  to  M  if,  for  every  ?  z  M  and  y  e  R  ,  <)>(y+0  31 
(y) .  An  important  invariance  property  of  distributions  derived  from 
Sjj(0,Ijj)  is  the  following;  a  proof  is  given  in  Jensen  and  Good  (1981). 

M 

LEMMA  2.  Let  L(Y)  be  spherically  symmetric  about  9  z  M,  and  let  R*  -*• 

R  be  location- invariant  with  respect  to  M  and  be  homogeneous  of  degree 
zero.  Then  the  distribution  L(4 >(Y))  on  R  is  invariant  for  all  distribu¬ 
tions  L(Y)  in  the  class  (S  (9,1  );  9  z  M) . 

N  -  -  N  '* 

Let  St  be  countable  and  let  (U(t);  t  e  SI}  be  a  stochastic  process  with 
values  in  R*3  such  that,  for  each  k  and  {t,  ,t_ . tt  }  z  M,  the  joint  dis- 

i  4.  iC 

tribution  of  {U(t^) , . . . ,U(t^) }  considered  as  a  distribution  on  R^  is 
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spherical;  then  {U(t);  t  e  W}  is  called  a  spherical  process .  A  count¬ 
able  process  is  spherical  if  and  only  if  it  is  a  scale  mixture  of  spherical 
Gaussian  processes,  i.e.  if  for  each  k  and  {t^.t^, • • • »t^}  e  M,  the  joint 

probability  density  function  (pdf)  of  U  =  [U(tp . U(t^)]  of  order  (p*k) 

is  given  by 

f (U)  =  (2r)"pk/2/^r~pkexp(-tr  0 ’U/2r2)dH(r)  (2.1) 


with  H(*)  a  mixing  distribution  on  (0,°°) ;  cf.  Hartman  and  Wintner  (1940). 
Various  choices  for  H(*)  yield  different  spherical  processes,  including 

the  spherical  Gaussian  process. 

2 

In  the  context  of  T  charts,  we  assume  that  successive  vector  obser¬ 
vations  are  generated  throughout  the  monitoring  period  from  a  spherical 
process  and  thus  are  dependent  except  in  the  Gaussian  case.  In  particular, 
partition  M  into  disjoint  subsets  corresponding  to  successive 

sampling  periods  with  a  base  period  of  length  m  and  succeeding  periods 
of  length  n.  Let  y(t)  and  £(t)  be  step  functions  on  M,  constant  within 
samples  with  possible  jumps  between  samples,  such  that  u(t)  e  RP  and 
Z(t)  (pxp)  is  positive  definite.  Our  model  for  generating  observations 
from  a  spherical  process  (U(t);  t  e  V)  with  values  in  Rp  is 


YQj  -  S!s(tj)U(tj)  +  v(tj);  t.  e  NQ,  j  -  1,2, ...m; 
=  ris(tj)U(tj)  +  |j(tj);  tj  e  W±,  j  *  1,2,.  ...n; 


(2.2) 

(2.3) 


where  i  =*  1,2,...  .  In  monitoring  location  parameters  the  usual  assump¬ 
tion  is  that  scale  parameters  remain  stationary,  i.e.  that  I(t)  =  Z  for 
every  t  e  N.  We  return  to  this  point  later. 


J 
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2 

3.  Properties  of  T  Charts.  In  this  section  we  study  Hotelling's 
2 

(1947)  T  charts  under  the  standard  validating  assumptions.  Specifically, 

we  establish  the  propensity  of  these  charts  to  signal  when  a  process  is  not 

in  control,  using  stochastic  comparisons  of  run  lengths. 

Let  Y.  »  fY. ,  ,Y. ,Y.  ]  be  the  vector  values  observed  during  the 
~i  ~il  ~i2  ~m 

ith  sampling  period,  and  let  TCY^)  be  the  statistic  to  be  charted  using  a 
procedure  for  monitoring  the  means  of  the  sampled  process.  In  testing 

2 

H:  y  =  Uq  against  A:  y  #  Pq  using  a  single  sample,  it  is  known  that  the  T 
test  (i)  is  invariant  under  nonsingular  linear  transformations  of  the  data 
and  thus  is  free  of  the  arbitrary  coordinates  used  to  represent  the  obser¬ 
vations,  and  (ii)  depends  on  Y^  only  through  Y^  and  S^,  and  its  power  de¬ 
pends  on  the  parameters  only  through  the  noncentrality  parameter  \  - 
n(y-yg)'2  '■(u-Pq).  Henceforth  let  x^  be  the  class  of  all  procedures  moni¬ 
toring  at  level  a  using  statistics  (T(Y^) ,T(Y2) , . . . }  having  the  property 
(i),  and  let  x2  be  the  class  of  procedures  monitoring  at  level  a  on  each 
occasion  and  having  property  (ii) .  In  the  case  of  a  single  sample  these 

classes  provide  a-level  tests  for  H:  y  =  y^  against  A:  y  ^  y^,  and  it  is 

2 

widely  known  that  the  T  test  is  uniformly  most  powerful  among  tests  in 

the  classes  x^  and  x2  on  a  fixed  sampling  occasion  (cf.  Anderson  (1958), 

Section  5.5).  The  notions  of  level  and  power  do  not  carry  over  to  the 

entire  monitoring  period,  however,  as  control  charts  typically  signal  with 

unit  probability  even  when  the  process  is  in  control.  Nonetheless,  we  next 

2 

establish  that  the  T  chart  is  optimal  among  procedures  in  the  classes  x^ 

and  x2  in  the  strong  sense  that  its  run  length  is  stochastically  smallest 

2 

at  a  fixed  alternative  y  ^  Uq.  The  T  chart  thus  tends  to  signal  more 
quickly  than  other  procedures  in  these  classes  when  the  process  is  not  in 
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control. 

THEOREM  L.  In  monitoring  the  means  of  a  p-dimensional  Gaussian  process, 

suppose  the  process  is  stationary  with  means  p  #  Uq  and  dispersion  matrix 

E.  Of  all  invariant  procedures  monitoring  at  level  a  in  the  class  t^,  and 

of  all  a-level  monitoring  procedures  in  the  class  t2»  the  run  length  of  the 
2 

T  chart  is  stochastically  smallest. 

Proof.  By  stationarity,  the  mutual  independence  of  {T(Yp ,T(Y  ) , . . . }  and 

their  identical  distributions,  it  follows  for  all  procedures  in  and 

that  their  run  lengths  are  geometric  with  parameter  equal  to  the  probability 

2 

of  exceeding  a  control  limit  on  any  particular  occasion.  Because  the  T 
test  is  optimal  among  a-level  tests  in  and  t 2>  the  probability  of  ex¬ 
ceeding  che  respective  control  limit  at  a  fixed  y  /  is  at  least  as  large 
2 

for  T  as  for  any  other  procedure  in  nr  The  proof  is  completed  on 

recalling  that  the  family  {G(t;8>;  8  s  (0,1)}  of  geometric  distributions 

2 

decreases  stochastically  in  8,  thus  assuring  that  the  run  length  of  the  T 
chart  is  stochastically  minimal. 

We  turn  next  to  modifications  of  the  standard  charts  which  might  be 
made  in  response  to  the  practical  difficulties  of  implementing  those  pro¬ 
cedures. 

2  2 

4.  Modified  Charts.  Both  the  standard  X  and  T  charts  require  that 

.  2 
the  control  value  Uq  be  known,  and  in  addition  the  X  chart  requires  the 

value  of  Z.  Often  these  parameters  are  unknown.  In  the  one-dimensional 

case  it  is  common  practice  to  estimate  such  parameters  in  a  base  period 

when  the  process  is  in  control,  and  to  use  these  estimates  in  lieu  of  the 

unknown  parameters.  Carried  over  to  the  multidimensional  case,  this  approach 
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yields  modifications  to  the  standard  charts  as  follows.  Let  Yg  =  ^Yoi,Y02’ 
. . . ,Y_  1  e  F  consist  of  sample  observations  in  a  base  period  of  length 

-Om  pxm  to 

m  yielding  one  or  both  the  estimates,  Yg  for  Uq  and  Sg  for  Z.  When  Z  is 

unknown  it  is  assumed  that  m  >  p  to  ensure  the  invertibility  of  Sg.  The 

modifications  of  the  standard  charts  to  be  considered  here  are  set  forth  in 

Table  1,  where  the  listed  statistics  are  appropriate  versions  of  (1.1)  and 

(1.2)  together  with  updated  versions  of  their  control  limits.  The  latter 

are  chosen  to  assure  that  the  modified  charts  will  continue  to  monitor  at 

level  a  on  each  sampling  occasion. 

Note  that  the  Type  2  and  Type  3  modifications  may  be  considered  as 
2  2 

substitutes  for  both  X  and  T  in  (1.1)  and  (1.2).  Practical  advantages  of 

2 

the  Type  2  chart  over  the  standard  T  chart  are  (i)  only  one  matrix  inver¬ 
sion  is  needed  using  the  Type  2  modification,  and  (ii)  sample  sizes  after 
the  base  period  need  not  exceed  p.  This  of  course  assumes  stationarity  of 
the  process  dispersion  parameters,  which  may  be  checked  periodically  using 
the  procedures  of  Hotelling  (1947)  or  of  Montgomery  and  Wadsworth  (1972). 

Run-length  distributions  of  these  modified  charts  generally  are  in¬ 
tractable;  they  are  not  geometric  owing  to  dependencies  among  the  success¬ 
ive  statistics  charted.  However,  the  main  result  of  this  section  is  that 
all  these  distributions  may  be  bounded  below  in  terms  of  standard  geometric 
distributions.  That  the  modified  charts  eventually  do  signal  is  shown  in 
the  proof  of  the  theorem. 

THEOREM  2.  Denote  the  typical  modified  statistic  of  Table  1  by  T(Y^,Yg) 
and  its  upper  control  limit  by  c^;  suppose  {Y^.Y^,...}  are  observations 
from  a  stationary  process;  and  let  N  be  the  run  length  of  the  modified 
chart.  Then  for  every  positive  integer  t, 
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P(N>t)  >_  1  -  G(t;B) 

where  3  is  Che  marginal  probability  3  =  PfT^.Yg)  >  c  ]  . 

Proof.  Let  M  be  geometric  with  distribution  G(t;3);  fix  t;  and  for  this 
value  note  that  P(M>t)  =  (1-8) C.  Similarly  express  P(N>t)  in  terms  of  the 
finite-dimensional  distribution  of  (T(Y^ ,Yq) .... ,T(J  ,Yq) } ,  and  apply  Lemma 
1  to  infer  that 

P(N>t)  =  P[T(Y  Y  )<c  ,...,T(Y  .Yj<c  ] 

t  U  a  ~t  ~u  a  (4.  ].) 

>  i21P[T(Yi,Y0)<Ca]  =  (1-3)  =  P(M>t) . 

The  proof  using  these  finite-dimensional  results  will  be  complete  if  it  can 
be  shown  that  the  modified  charts  eventually  signal  with  unit  probability. 

To  see  this,  use  the  expression  given  for  P(M>t)  and  the  conditional  inde¬ 
pendence  of  (TCY^Yq)  , . . .  ,T(Yt,Y0)}  given  YQ  to  write 

P(N>t)  =  / F  (P [T(Y,YQ) <0^ | Yq] } CdF(Y0)  (4.2) 

pxm  'J 

where  F(*)  is  the  cdf  of  Yn  and  Y  has  the  same  distribution  as  Y.  for  i  = 

~0  -i 

l,2,...,t.  Clearly  the  function 

H(t,YQ)  -  {P[T(Y,Y0)<Ca|Y0]}t  (4.3) 

satisfies  H(t,Yg)  <  1  and  H(t,Yg)  ■*  0  as  t  -*•  «  for  each  fixed  Y^.  Moreover, 

/  dF(Yg)  =  1.  From  the  dominated  convergence  theorem  it  follows  that 
'p><m 

£im  P(N>t)  -  /F  £im  H(t,Y0)dF(Y0)  -  0.  (4.4) 

pxm 

This  shows  that  the  modified  charts  eventually  signal  with  unit  probability 
and  thus  that  the  proof  holds  for  every  t  >  0. 


We  note  that  Theorem  2  holds  under  weakened  assumptions.  The  role  of 
Gaussian  assumptions  in  that  theorem  is  to  supply  the  control  limit  c^;  the 
theorem  holds  otherwise  as  long  as  {Yq,Y^ . . . }  are  mutually  independent 
and  {T(Y^ ,Yq) .TCY^jYq) , . . . }  are  identically  distributed.  In  particular,  sta- 
tionarity  is  not  required  for  the  Type  4  modification  as  long  as  the  non¬ 
centrality  parameters  {X^  =  [nm/ (n+m) K^-Uq)  ' ^ (y^-Ug)  ’  i  =  are 

held  constant. 

5.  Drifting  Processes.  Thus  far  we  have  supposed  that  the  monitored 

process  is  stationary  beyond  a  base  period,  whether  or  not  it  is  in  control. 

Combined  with  independence,  stationarity  leads  to  charts  whose  run  lengths 

2  2 

are  geometric.  This  section  develops  properties  of  the  X  and  T  charts  in 

monitoring  the  means  of  nonstationary  multidimensional  processes  which  we 

call  drifting  processes.  The  run  lengths  for  drifting  processes  are  not 

geometric  even  under  independence  because  the  successive  statistics  are  not 

distributed  identically.  Such  distributions  are  complicated  further  by  their 

dependence  on  a  countable  sequence  of  parameters.  Fortunately,  intractable 

distributions  of  this  type  can  be  bounded  by  simpler  distributions  under 

suitable  conditions  on  the  sequence  of  parameters. 

Let  5  =  (6^,62....)  be  a  sequence  of  parameters  for  a  drifting  process; 

let  F ( t ; 5 )  be  the  distribution  of  run  lengths  for  the  corresponding  chart; 

and  let  V  be  the  set  of  all  bounded  sequences.  Two  sequences  §  and  5*  in 

V  are  said  to  be  ordered  as  5  <_  5*  whenever  6^  <_  5*  for  all  i.  Under  this 

partial  ordering  an  element  6  e  P.  is  said  to  be  minimal  for  P.  e  P,  and 

~m  U  -  U 

5W  t  Vn  to  be  maximal  for  P„,  if  6  <  5  <  5„  for  every  5  t  Vn.  The  family 

0  -  0  -m  —  ~  —  ~M  -  0 

{F(t;5);  S  z  V}  of  distributions  is  said  to  be  stochastically  decreasing 
in  6  if,  for  any  5  and  5*  in  V,  the  ordering  5  <_  5*  implies  the  stochastic 


ordering  F(t;5*)  _>  F(t;6)  for  every  t.  For  drifting  processes  this  prop¬ 
erty  captures  the  notion  that  a  chart  will  signal  more  quickly  as  the  pro¬ 
cess  drifts  further  from  control. 

An  important  property  of  a  stochastically  decreasing  family  of  run- 

length  distributions  is  that  an  envelope  of  curves  can  be  constructed  for 

certain  members  of  the  family  in  terms  of  simpler  distributions.  Given  a 

bounded  sequence  8  =  (5,,S„,...),  let  8  =  inf  ($,  .)  and  ■  sup(<5,  , 

~  1  Z  m  L  Z  ”  i 

5,,...),  and  define  o(m)  =  (5  ,6  ,...)  and  6(M)  =  (6W»5W, ...).  A  basic 
Z  ~  mm  ~  MM 

result  is  the  following. 

LEMMA  3.  Suppose  the  family  (F(t;5);  5  e  V]  of  run-length  distributions  is 
stochastically  decreasing  in  6.  Then  for  any  8  t  V,  the  inequalities 

F(t;6(M))  >  F(t;6)  >  F(t;fi(m)) 

hold  for  every  t  >  0.  Moreover,  these  bounds  hold  for  every  F(t;6*)  such 
that  S(m)  <_  6*  <_  5(M). 

Proof.  The  first  conclusion  is  a  consequence  of  the  transparent  ordering 
5(m)  <_  6  <_  <5(M)  and  the  fact  that  the  family  {F(t;5);  5  e  V}  is  stochasti¬ 
cally  decreasing  in  6.  The  second  conclusion  follows  on  repeating  these 
arguments  for  any  6*  in  <5(m)  <_  5*  5(M). 

It  may  be  noted  that  F(t;6)  depends  on  the  path  5  = 
the  drifting  process,  whereas  the  upper  and  lower  bounds  do  not.  This 
fact  simplifies  a  study  of  the  bounding  distributions,  which  are  geometric 
in  some  important  cases. 

Lemma  3  supports  bounds  on  distributions  of  run  lengths  for  the  stan- 
2  2 

dard  X  and  T  charts  and  their  modifications  under  drifting.  What  is 
needed  in  each  case  is  to  identify  the  appropriate  parameter  sequence  8 
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and  to  show  that  {F(t;5);  5  e  V]  is  stochastically  decreasing  in  5. 

2 

To  fix  ideas,  first  consider  the  X  chart  with  run  length  N  under  a 
drifting  process  with  parameters  {(y^.Ej),  » *  *  •  f°H°ws  that 


P(N>t)  =  ^PCX^cJ 


(5.1) 


2  2 

where  L(X.)  =  x  (p,^),  the  noncentral  chi-squared  distribution  with  p 

degrees  of  freedom  and  the  noncentrality  parameter  6  =  n^^-iO 'E^Cu^-Uq)  . 

The  parameter  sequence  for  F(t;6)  accordingly  is  6  =  (d^,^****)*  an<^  be- 
2 

cause  PCX^c^)  is  a  decreasing  function  of  5^,  it  follows  that  (i)  the 

family  {F(t;5);  5  e  V}  is  stochastically  decreasing  in  5,  and  (ii)  F(t; 

o (m) )  =  G(t;  a(m) )  and  F(t;  6(M))  =  G(t;  a(M)),  where  ct(m)  and  a(M)  are  the 

probabilities  of  exceeding  the  control  limits  on  those  occasions  for  which 

5.  *  5  and  o.  «*  5„,  respectively.  Similar  developments  apply  to  the  stan- 
i  m  j  M 

2  2 
dard  T  chart  using  the  monotonicity  of  the  noncentral  T  distribution  in 

its  noncentrality  parameter  5^,  leading  to  the  following. 


THEOREM  3.  Let  F(t;5)  with  5  e  V  be  the  run-length  distribution  of  either 
2  2 

the  standard  X  or  T  chart  under  a  drifting  process  with  parameter  sequence 
{5i  *  hCu^Uq) 'X~1(ui-U0) ;  i  *  1,2,...}.  Then  the  stochastic  bounds 

G(t;o(M))  >  F(t;5)  >_  G(t;a(m)) 


apply  for  every  £  >  0,  where  a(m)  and  a(M)  are  probabilites  of  exceeding 
the  control  limits  when  <5^  -  Sm  =  inf  . . .)  and  6^  =  6^  *  supCd^.d^, 
...),  respectively. 


Further  conclusions  are  evident.  Because  0  =  (0,0,...)  e  V  and  is 
minimal  for  V,  it  follows  that  the  geometric  distribution  G(t;a)  arising 
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when  the  process  is  in  control  stochastically  dominates  the  run-length 

distribution  for  any  drifting  Gaussian  process.  Because  the  former  signals 

2  2 

with  unit  probability,  it  follows  that  the  standard  X  and  T  charts  even¬ 
tually  signal  under  drifting.  The  behavior  of  the  standard  charts  may  be 
compared  under  two  drifting  processes  in  terms  of  5^  *  nOa^-Ug)  . 

When  =  Z  for  all  i,  the  inequality  6^  >_  5*  =  n(u£-Pg)£  ^  (ti u o>  asserts 

that  u.  is  more  distant  from  u_  than  u*  in  a  non-Euclidean  metric.  The 
-l  ~u  ci 

stochastic  ordering  F(t;6)  F(t;S*)  assures  that  the  chart  tends  to  signal 

more  quickly  as  the  process  drifts  further  from  control.  Even  if  a  process 
is  stationary  in  its  means,  i.e.  =  u  for  all  i,  the  charts  benefit  from 
successively  tightened  dispersion  parameters  through  refinement  of  the  pro¬ 
cess.  Specifically,  let  and  be  the  disPersion 

matrices  of  two  processes  such  that  Z*  -  is  positive  semidefinite  for  all 
i.  It  follows  that  5^^  >_  <5*,  i.e. 

n (y— Uq)  ’I^(ij-Uq)  >_  n(u-y0) '(E*)  Vu-Uq).  (5.2) 

hence  that  F ( t ; 5 )  F(t;5*)  for  every  t  >  0,  and  thus  that  the  chart  for  the 

former  process  signals  more  frequently  than  the  latter  at  a  fixed  u  #  Uq. 

Corresponding  results  hold  for  modified  charts  as  developed  in  Section 
4.  If  N  is  the  run  length  of  a  modified  chart  we  have,  in  the  notation  of 
Section  4, 

P(N>t)  -  P(T(Y1,Y0)<ca,...,T(Yt,Y0)<Ca) 

•  h  JlHT^.YpicjY  JdFOT  > 

pxm 

for  every  positive  integer  t.  For  modified  charts  of  Type  1  and  Type  4  the 
family  (F(t;5);  o  z  V)  of  run-length  distributions  is  stochastically 
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2  2 

decreasing  in  6,  conditionally  using  properties  of  X  and  T  charts  as  in 
developments  leading  to  Theorem  3,  and  unconditionally  because  the  in¬ 
equality  holds  point-wise  at  each  fixed  Y^.  For  Type  2  and  Type  3  modifi¬ 
cations  the  needed  ordering  is  supplied  in  the  following  lemma  under  the 
assumption  that  the  process  dispersion  parameters  are  stationary.  Recall 
that  these  procedures  otherwise  would  not  be  used. 


LEMMA  4.  Suppose  the  process  dispersion  parameters  are  stationary.  Then 

the  family  (F(t;6);  <5  e  V}  of  run-length  distributions  is  stochastically 

2 

decreasing  in  5  for  both  the  modified  T  charts  of  Type  2  and  Type  3. 

Proof.  For  Type  3  charts  temporarily  fix  Y^  and  note  that  these  essentially 
are  of  Type  2  conditionally.  For  Type  2  charts  we  have 


P(N>t)  *  P(T(Y1,Y0)<Ca,...,T(Yt,Y0)<Ca)  (5.4) 

where  T(Y^,Yq) ,  as  defined  in  Table  1,  may  be  rewritten  as 

T(Yi,Y0)  =  na^rs"1^-^)/  w.  (5.5) 

with 


"i  '  <L-“o,t'l(L-to)/<?i-So)'§ol(L-So> 


(5.6) 


2 

for  i  »  l,2,...,t.  A  standard  result  is  that  (.(w^  =  x  (m-p+1)  indepen¬ 
dently  of  Y^.  Write  w  *•  (w^ .w^, . . . .w^)  and  let  G(w)  be  their  joint  cdf. 
Now  recalling  the  definition  of  X^,  writing  (5.5)  as  TCY^.Yq)  *•  X^/w^,  and 
noting  that  {T(Y^ ,Yq) , . . . ,T(Yt ,Yq) }  are  conditionally  independent  given 
w,  we  evaluate  (5.4)  as 


u  2 

P(N>t)  -  fgt  ^PU^WjcJ^dGtw) 


(5.7) 
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t  t  2 

where  R+  is  the  positive  orthant  of  R  .  Because  conditionally  has  a  non¬ 
central  chi-squared  distribution  with  noncentrality  parameter  5^  *  nCy^-u^) ’ 

“1  2 
Z  ty^-y^),  the  conditional  probability  P(X^£w^ca|w)  decreases  as  <5^  in¬ 
creases,  point-wise  for  each  fixed  w^,  i  =  l,2,...,t.  For  Type  2  charts 
the  desired  property  holds  unconditionally  on  taking  expectations  using 
(5.7).  For  Type  3  charts  these  conditional  results  hold  point-wise  for  each 
fixed  Yq  and  thus  unconditionally  using  a  standard  argument. 

The  foregoing  results  now  may  be  combined  with  Lemma  3  to  give  the 
following. 

THEOREM  4.  Let  F(t;5)  with  &  e  t?  be  the  run-length  distribution  of  a  modi¬ 
fied  chart  under  a  drifting  process  with  parameters 

(i)  3  nCy^UQ) ’Ei1(ui-y0> ;  i  =  1,2,...}  for  Type  1  and  Type  4 
charts ,  and 

(ii)  {6i  -  n^-y^'E  ^v^-Uq);  i  =  1,2,...}  for  Type  2  and  Type  3 
charts.  Then  the  bounds 

F ( t ; 6 (M) )  >  F(t;6)  >  F(t;6(m)) 

hold  for  every  positive  integer  t,  where  5(m)  =  (6  ,  <5  ,...)  and  <5(M)  = 

~mm 

(<5jl»  •  •  • )  Mith  (5in  =  infCdj.d  ,...)  and  5M  =  sup(51 ,5^, . . . ) . 

It  is  of  interest  that  F(t;5(m>)  and  F(t;6(M))  are  precisely  distribu¬ 
tions  of  the  types  considered  in  Section  4.  Theorem  2  accordingly  yields 
the  bounds 


F(t;6(m))  <_  G(t;a(m))  and  F(t;6(M))  £G(t;a(M)) 


where  a(m)  and  a(M)  are  the  probabilities  of  exceeding  the  upper  control 
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limit  on  those  occasions  for  which  6.  =  6  and  5.  =  5W,  respectively. 

i  m  j  M 

It  was  noted  earlier  that  envelopes  for  the  run-length  distributions 

of  certain  drifting  processes  may  be  constructed  using  geometric  distribu- 

2 

tions.  To  illustrate,  envelopes  are  given  in  Figure  1  for  the  standard  X 
2 

and  T  charts,  each  monitoring  at  level  a  =  0.05  with  p  *  2  and  n  =  6.  The 
case  6=0  applies  when  the  monitored  process  is  in  control.  For  each  type 
of  chart  the  curves  labeled  6  =  3.0  and  5  =  6.75  represent  stationary  pro¬ 
cesses  having  the  parameter  sequences  5^  =  (3.0,3. 0, .. .)  and  6^  =  (6.75, 
6.75,...),  respectively.  These  curves  comprise  an  envelope  containing  the 
distributions  of  run  lengths  for  all  drifting  processes  satisfying 

3.0  <  6  <  6.  <  6U  <  6.75 

—  m  —  i  —  M  — 

2  2 

using  X  or  T  charts  as  appropriate.  Entries  in  Figure  1  were  found  on 
2  2 

converting  X  and  T  into  F-statistics  and  using  P.C.  Tang's  tables  as 

given  in  Graybill  (1961)  to  determine  the  probability  S  of  surpassing  the 

a-level  control  limit  on  any  particular  occasion.  Values  of  G(t;8)  then 

were  computed  directly  for  various  values  of  t  and  were  graphed. 

2  2 

It  is  informative  to  compare  curves  for  the  X  and  T  charts  at  a 

2 

fixed  u  #  Then  the  run  length  of  the  X  chart  is  stochastically  small- 

2 

er  than  that  of  the  T  chart.  This  property  may  be  observed  in  Figure  1 

for  the  cases  treated  there,  and  it  can  be  shown  analytically  for  any 

2 

choice  of  p,  n,  and  y  ^  y^.  The  greater  efficiency  of  using  X  when  appro¬ 
priate  on  any  monitoring  occasion  thus  translates  into  a  greater  propensity 

2 

to  signal  that  y  +  y q  in  comparison  with  the  T  chart. 

2 

The  optimality  of  the  standard  T  chart  among  procedures  in  and 
was  established  in  Section  3  for  stationary  processes  not  in  control.  This 
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optimality  carries  over  to  drifting  processes  as  well. 

THEOREM  5.  Consider  a  p-dimensional  Gaussian  process  with  drifting  para¬ 
meters  » •  •  •  ^  on  successive  sampling  occasions,  to  be 

monitored  for  means.  Of  all  invariant  procedures  monitoring  at  level  a  in 

the  class  x ^ ,  and  of  all  ct-level  monitoring  procedures  in  the  class  t^7  the 
2 

run  length  of  the  T  chart  is  stochastically  smallest. 

Proof.  Let  {T(Y,),  T(Y„),...}  be  the  statistics  and  c  the  control  limit 
-  - 1  ~2  a 

of  any  procedure  in  x^  or  and  let  N  be  its  run  length.  Because 

t 

P (N>t)  =  H  P(T(Y.)<c  )  , 

1*  J-  ■*'1  Ot 

it  is  clear  that  (i)  the  run-length  distribution  F(t;3)  of  the  chart  de¬ 
pends  on  the  parameters  §  =  (g^g,,,...)  with  {0i  =  PCTCY^c  );  i  =  1,2, 
and  the  family  {F(t ; 0) ;  g  e  $}  decreases  stochastically  in  §.  Let 

3  be  the  sequence  of  parameters  for  any  procedure  in  x^  or  x^  and  let  6* 

2 

be  the  corresponding  sequence  for  the  T  chart.  Because  on  each  sampling 
2 

occasion  T  is  optimal  among  procedures  in  x^  and  x^,  it  follows  that 
3  6*.  The  conclusion  of  the  theorem  now  follows  from  the  stochastically 

decreasing  character  of  {F(t;S) ;  g  e  V} . 

6.  Monitoring  Spherical  Processes.  Many  findings  of  the  foregoing 
sections  carry  over  to  observations  generated  from  a  spherical  process 
{U(t);  t  s  N}  as  in  Section  2.  Consider  the  standard  X  and  T  charts 
using  outcomes  {Y  (t)}  of  the  vector-valued  process  defined  in  (2.3) 
having  the  location  parameters  y(t)  and  the  dispersion  parameters  I(t). 

For  convenience  we  drop  the  arguments  and  index  these  parameters  as  (y^, 

5^)  for  the  ith  monitoring  period.  The  principal  properties  of  these 
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charts  under  stationary  and  drifting  processes  induced  from  a  spherical 
process  are  the  following. 


THEOREM  6.  Let  c  (X)  and  c  (T)  be  a-level  control  limits  for  the  X  and 
-  a  a 

2 

T  charts  for  use  with  observations  generated  from  a  spherical  process. 

Let  N  and  N  be  the  run  lengths  of  these  charts  having  the  distributions 
X  1 

Fx(t;5)  and  FT<t;6)  with  parameters  =  n(yi-yQ) ,E11 (ui-uQ) ;  i  =  1,2,...}, 


(i)  c  (T)  »  T  (p,n-l) ; 

a  a 

(ii)  P(NT<t)  =  G(t;a)  when  y(t)  *  for  all  t; 

(iii)  P(N  >t)  >  1  -  G(t;8)  when  the  process  is  stationary,  where 

A  1 

3  =  P(X2>c  (X));  and 
a 

(iv)  the  families  {Fx(t;<S);  {  £  5}  and  { FT ( t ; 5 ) ;  5  e  V}  are  stochasti¬ 
cally  decreasing  in  5  *  (^1*^2* **')* 

Proof,  (i)  Starting  with  Y  *  [Y  ^  ,Yi2’  ‘  *  ’  ,Yin^  ’  make  tlie  location-scale 

-U 

changes  Z,  =  [Z. . ,Z.0». . . ,Z.  ]  with  Z..  *  I.  (Y..-y„)  and  note  that 

-1  -ii  -iz  -in  ~ij  -I  -0 


(6.1) 


where  Z ^  and  respectively  are  the  sample  mean  vector  and  the  sample  dis¬ 
persion  matrix  computed  from  Z^.  Moreover,  (.(Z^)  is  spherical  on  F 

Thus  choosing  M  as  {0}  c  F  and  noting  that  4>(Z.)  ■  nZ.’W  ^Z.  is  homogen- 

-*  pxn  ~ l  '■'i**!  ~ l 

2 

ous  of  degree  zero  in  its  argument,  we  apply  Lemma  2  with  (Z^)  =  T^  to 
2  2 

infer  that  L(T^)  ■  T  (p,n-l)  for  any  underlying  spherical  process  when 

2 

y(t)  ■  un  for  all  t.  The  conclusion  c  (T)  »  T  (p,n-l)  follows  immediately. 

v  w  Cl  Ct 

(ii)  For  each  fixed  t  write 


P(N _>t)  -  P(T?<c  (T) . T2<c  (T) ) ; 

T  a  t  a 
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define  Z  *  [Z.,Z  ,...,Z  ]  e  (F  ) t  in  terms  of  the  standardized  variables 
■*  ■“  J»  2,  “*  t  p^n 

introduced  in  the  proof  of  (i) ;  and  note  that  L(Z)  is  spherical  on  (Fpxn)C. 
Accordingly  take  M  as  {0}  =  ^pxn^’  0^serve  that  i>  (Z)  =  (nZ’W^Z^ , . . . , 

nZ^W^Z^  is  homogeneous  of  degree  zero  in  its  argument;  and  apply  Lemma  2 

2  2  2  2 

once  more  with  <£(Z)  =  (T^, _ ,Tt)  to  conclude  that  L(T^#...fT  )  is  iden¬ 

tical  to  its  normal-theory  form  when  p(t)  =  p^.  Conclusion  (ii)  now  follows 
directly  from  (6.1)  and  (6.2)  and  results  known  for  the  Gaussian  case. 

(iii)  From  the  representation  (2.1)  for  the  underlying  spherical  pro¬ 
cess,  P(N^>t)  may  be  written  as 

P(N_>t)  «  P(X?<c  (X),...,xf<c  (X)) 

A  1  Ct  L  Ot 

„  t  ,  (6.3) 

«  f0  i21p<Xif.Ca(X)|r)dH(r) 

2  2 

because  {X^,...,Xt>  are  conditionally  independent  given  r.  Under  the  hypo¬ 
thesis  of  stationarity  (6.3)  becomes 

P(Nx>t)  =  /”[P(X2<ca(X) |r)]tdH(r)  (6.4) 

2  2  2 
where  X  has  the  typical  distribution  of  (X^,...,X  }.  Conclusion  (iii)  now 

follows  from  the  arguments  leading  to  Lemma  1. 

(iv)  Consider  { F^( t ; 6 ) ;  <5  e  V} .  From  (6.3)  it  is  clear  that  the  con- 

—2  2  2 

ditional  distribution  L(r  X^|r)  is  x  (p,5^(r))  with  the  noncentrality 
parameter 

6t(r)  »  n(ui-uQ)  'z”1  (^-Uq) /r2.  (6.5) 

That  the  conditional  run-length  distributions  {F  ( t ; 6 | r) ;  5  £  V}  are  sto- 

A  'w 

chastically  decreasing  in  <5  follows  as  in  the  proof  of  Theorem  3  point-wise 
for  each  fixed  r.  The  property  holds  unconditionally  on  applying  a 
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standard  argument.  The  case  of  {FT(t;<S);  6  e  V)  is  treated  similarly. 

2 

Conclusion  (i)  assures  that  the  standard  T  chart  with  the  normal- 
theory  control  limit  c^  is  appropriate  for  monitoring  any  spherical  pro¬ 
cess  after  location  and  scale  changes,  at  level  a.  However,  the  value  c^(X) 

2 

required  for  monitoring  at  level  ot  using  the  X  chart  generally  depends  on 

the  underlying  spherical  process.  Conclusion  (ii)  of  Theorem  6  establishes 

2 

the  familiar  normal-theory  run-length  distribution  of  the  T  chart  for  any 
spherical  process  when  the  process  is  in  control. 

Results  for  modified  charts  given  in  Sections  4  and  5  under  Gaussian 
assumptions  carry  over  to  spherical  processes.  Proofs  for  these  extensions 
follow  the  pattern  of  proof  for  conclusions  (iii)  and  (iv)  of  Theorem  6. 

In  view  of  the  representation  (2.1),  conditional  properties  already  have 
been  established  in  the  earlier  sections  point-wise  for  each  r;  uncondition¬ 
al  properties  then  yield  to  a  standard  argument.  Further  details  are  omit¬ 
ted  in  the  interest  of  brevity. 

Our  final  results  deal  with  the  relative  performance  of  a  given  chart 

under  alternative  spherical  processes.  By  a  given  chart  is  meant  either 
2  2 

the  X  or  T  chart  with  a  fixed  value  for  its  control  limit.  Let  p(*)  and 
v(*)  be  finite-dimensional  measures  characterizing  two  zero-mean  stochastic 
processes.  If  for  every  N,  u(')  is  more  peaked  about  0  e  R  than  v(*)  in 
the  sense  of  Sherman  (1955),  then  the  p-process  is  said  to  be  more  peaked 
about  the  zero  function  than  the  v-process.  In  the  case  of  spherical  pro¬ 
cesses  having  the  representation  (2.1),  it  can  be  shown  that  of  two  such 
processes  having  the  mixing  distributions  H  (r)  and  H^(r),  the  p-process 

is  more  peaked  about  zero  than  the  v-orocess  if  and  only  if  H  (r)  >  H  (r) 

p  —  v 

for  each  r  >_  0.  Our  principal  result  is  the  following. 


I 
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THEOREM  7.  Lee  F  (t;5)  and  F  (t;<5)  be  the  run-length  distributions  of  a 

-  H  -  \) 

chart  under  drifting  processes  having  the  same  parameters 

i5i  =  n(yi_!fo) '?il(Hi“y0) ;  1  = 

and  generated  from  spherical  processes  having  the  measures  y(*)  and  v(*)> 
respectively.  If  the  u-process  is  more  peaked  about  zero  than  the  v-process, 
then 

2 

(i)  F^(t;0)  >_  F^(t;0)  for  the  X  chart  when  the  processes  are  in  con¬ 
trol;  and 

2 

(ii)  F  ( t ; 5 )  >  F  (t;£)  for  the  T  chart. 

u  -  —  v 

2 

Proof.  (i)  Let  c(X)  be  the  control  limit  for  the  X  chart,  and  let 

Ai  =  {?ieFpxn!Xi=n?i?ilc(X)}  (6-6) 

in  the  notation  used  in  the  proof  of  Theorem  6.  On  letting  A  =  A^nA^. . .nAt, 
we  infer  from  the  first  expression  on  the  right  of  (6.3)  that 

[P(Nu>t)-P(Nv>t)]  =  [jj  (A)-v  (A)  ]  .  (6.7) 

But  A  is  a  convex  subset  of  (F  )C  symmetric  about  0  under  reflection. 

pxn 

Under  the  hypothesis  that  u(*)  is  more  peaked  about  0  than  v(-),  it  follows 

that  [p (A)-v (A) ]  >_  0,  which  is  equivalent  to  conclusion  (i)  . 

2 

(ii)  Consider  the  standard  T  chart  with  control  limit  c  .  Let  X  be 

a 

a  typical  run  length  and  let  and  be  the  run  lengths  of  the  chart 
under  the  two  processes.  For  any  underlying  spherical  process  note  as  in 
(6.3)  that 

P  (N>t)  =  P(T^±ca . T^<ca) 

-  /q  iHlP(T^<ca!r)dH(r) 


(6.8) 
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2  2  2 

where  {T^ . . .  ,Tt)  are  conditionally  independent  given  r,  and  the  con- 

2 

ditional  distribution  KT^jr)  is  a  noncentral  Hotelling's  (1931)  distri¬ 
bution  depending  on  r  only  through  the  noncentrality  parameter  ^(r)  given 
by  (6.5).  Define 


G(t;S(r)) 


iV^Jr) 


(6.9) 


2 

with  <5 (r)  =  (5^ ,6^, . . . ,<5t) /r  .  Applying  the  relation  (6.8)  twice  and  using 

(6.9),  we  evaluate  the  difference  D(t)  =  [P(N  >t)-P(N  >t)l  as 

v  U 


D(t)  =  /^G(t;d(r))dHv(r)  -  /^G(t;5(r)dH^(r) . 


(6.10) 


Under  the  hypothesis  that  H^(*)  dominates  H^(0  stochastically,  there  are 
increasing  functions  ?(v)  and  n(v)  with  C(v)  £  n(v)  and  a  random  variable 
V  with  distribution  H(-)  such  that  JL($(V))  -  H^(*)  and  L(ri(V))  =  H^O); 
cf.  Lehmann  (1959),  page  73,  It  follows  that 


D(t)  =  j^[G(t;5(n(v)))  -  G(t;o(£(v)))]dH(v) .  (6.11) 

But  G(t;6(r))  is  an  increasing  function  of  r;  because  n(v)  >_  £(v)  for  each 
v,  it  follows  that  the  integrand  on  the  right  of  (6.11)  is  point-wise  non¬ 
negative  and  thus  the  integral  is  nonnegative.  This  gives  P(Ny>t)  >_P(N^>t), 
which  is  equivalent  to  conclusion  (ii)  of  the  theorem. 


It  was  noted  earlier  that  the  control  limit  ^(X)  for  monitoring  at 
2 

level  a  using  the  X  chart  generally  depends  on  the  underlying  spherical 

process.  Suppose  the  control  limit  c*(X)  is  chosen  for  some  reference  dis¬ 
ci 

tribution  such  as  the  Gaussian  distribution.  Provided  the  processes  are 
in  control,  the  run-length  for  the  reference  distribution  as  a  bound 
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stochastically  dominates  the  run  lengths  for  all  spherical  processes  less 

peaked  than  the  reference,  and  is  dominated  stochastically  by  the  run  lengths 

for  all  spherical  processes  more  peaked  than  the  reference  distribution. 

These  results  are  consequences  of  conclusion  (i)  of  Theorem  7.  Conclusion 

2 

(ii)  of  that  theorem  assures  that  the  T  chart  will  tend  to  signal  more 
frequently  under  drifting,  the  more  peaked  the  measure  of  the  underlying 
process. 
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TABLE  1.  Statistics  and  their  control  limits  in  modified  X  and  T  charts 
for  monitoring  the  means  of  a  multidimensional  production  process. 


Modification 

Type 

Statistic 

Control 

Limit 

1 

[tun/  (n-hn)  ]  (Y^Yq)  *  Z_1  (Y^-Yg) 

x>> 

2 

T^(p?m-1) 

3 

[nm/  (n-hn)  ]  (Y^Yq)  ’  S^1  (Yj-Yq) 

T^(pfm-1) 

4 

[ran/ (n-Hn)  ]  (Y.-Yg)  * s"1  (Yj-Yg) 

T^(p,n-1) 
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